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We argue that a multiband superconductor with sign-changing gaps may have multiple spin resonances. We 
calculate the RPA-based spin resonance spectra of a pnictide superconductor by using the five band tight-binding 
model or angle-resolved photoemission spectroscopy (ARPES) Fermi surface (FS) and experimental values of 
superconducting (SC) gaps. The resonance spectra split in both energy and momenta due to the effects of multi- 
band and multiple gaps in s*— pairing; the higher energy peak appears around the commensurate momenta due 
to scattering between a— FS to 7/ 5— FS pockets. The second resonance is incommensurate coming from /3— FS 
to 7/(5— FS scatterings and its vector is doping-dependent and hence on the FS topology. Energies of both 
resonances 10].^^ are strongly doping dependent and are proportional to the gap amplitudes at the contributing 
FSs. We also discuss the evolution of the spin excitation spectra with various other possible gap symmetries, 
which may be relevant when either both the electron pockets or both the hole pockets completely disappear. 

PACS numbers: 74.70.Xa,74.40.-n,74.20.Rp,74.25.Jb 



Introduction: The magnetic resonance behavior which is 
directly probed by the inelastic neutron scattering (INS) spec- 
troscopy gives valuable information about the pairing mech- 
anism of the unconventional superconductors. In cuprates, 
INS exhibits a clear signature of a resonance mode in addition 
to its characteristic dispersion (known as 'hour-glass' behav- 
ior) which is enhanced dramatically below Tc, and the mode 
energy scales universally with the SC gap amplitude. fTI In 
a multiband unconventional superconductor, the situation be- 
comes more complex due to the presence of the multiband and 
multiple SC gaps as well as the possibility of having multiple 
pairing symmetries. Given the important role the spin reso- 
nance played in the identification of pairing symmetry in the 
past, we investigate the details of spin resonance in a multi- 
band model relevant for pnictide superconductors. The pri- 
mary question we ask: Could multiple superconducting gaps 
at various bands lead to multiple spin resonances as opposed 
to the case of single spin resonance for single band super- 
conductors? We find that indeed multiband superconductors 
could have multiple spin resonances. It is also important to 
study how the resonance spectra evolve in (q, —phase space 
as function of doping and other parameters. 

An example of a multiband superconductor is MgB2. It is 
a conventional s— wave superconductor and its signal in the 
INS measurement is absent as the magnetic spectra are sen- 
sitive to the sign change of the SC gap on the FS sheet. fJl 
Recently discovered pnictides offer a new testing ground for 
these ideas. |[3]|4| In pnictides, it is now generally accepted that 
an unconventional pairing symmetry is present with evi- 
dence coming from ARPES [SJ and other probes [6, 7 1. This 
conclusion is also supported by theoretical calculations|l8l|9]. 

We have performed a standard BCS susceptibility calcula- 
tion within the RPA framework with s*— wave pairing sym- 
metry for the materials with the Fe-122 band structure. The 
computed dispersion and intensity of the spin polarized INS 
spectra reveal the following: (i) Two SC gaps set the scale 



for the two spin resonances lu},^ ,. and w^g^, that arise from the 
fact that the SC gap magnitudes (although the same pairing 
symmetries) are different at two hole pockets at F. At the 
same time the momenta for spin resonances are different: w^^^ 
is commensurate and w^^^ is incommensurate (with doping- 
dependent g— vectors), (ii) We also find that both w^g^ and 
w^gj, are also doping dependent but universally proportional 
to the SC gap amplitudes. When one of the hole pockets 
disappears, one of the resonances also disappears, (iii) Fur- 
thermore, we also show that when both hole pockets or both 
electron pockets vanish completely from the FS, the pairing 
symmetry should be changed to one of the form which will 
give a change of sign at the same FS and as a result the intra- 
band scattering spectra will be shifted to the low-energy and 
low -momenta scale. lID 

The large energy scale of both modes oj^^^ ~ 8 — 15meV 
and w^gg ~ 5 — 12meV depending on doping makes it likely 
that both resonances can be seen in INS. There is clear evi- 
dence for the first commensurate resonance lfTOl [TTI and there 
is also some evidence for the second mode lfTom2l whose en- 
ergy and q— values agree well with our calculation. 

Formalism: The bare dispersion and its corresponding or- 
bital characters are calculated by solving the five d— orbitals 
(in one Fe unit cell) tight-binding formalism. lfT3l We obtain 
doping by the rigid band shift method. In a multiorbital SC 
state, the BCS free-fermion susceptibility is a tensor[8J given 
by 

Xorst«(q, jw,„) = ^ M;;';„(k,q)[G'k(icj„) 

Here v {v') is the initial (final) eigenstate coming from s,t 
(r, u) orbitals and G and F are the normal and anomalous 
part of the Green's function, respectively. M is the orbital 
to band matrix element (ME) made of the eigenvectors. The 
many-body poles in the dressed x are obtained within the RPA 
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framework. The terms in the spin RPA interaction vertex Us 
that are included in the present calculation are intraorbital in- 
teraction U, an interorbital interaction [/', Hund's rule cou- 
pling J > 0, and the pair hopping strength J'. The rest of the 
interaction can be expressed in term of U U' = U ~ J/2, 
J — U/A and J' = J |8| While the calculations are done 
in the unfolded one Fe per unit cell Brillouin zone (BZ), the 
results are subsequently folded to the two Fe per unit cell BZ. 
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FIG. 1: (Color online) (a)-(b) RPA-BCS x"('^). plotted for calcu- 
lation including the ME effect, M from Eq. 1 in (a) and without 
including the ME effect {M — 1) in (b). (d)-(g) The correspond- 
ing lower panel figures clarify the mechanism of double resonances 
at two representative momenta cuts. x"s are shifted by a constant 
value to ease comparison with others, while they all go to zero at 
oj = 0. In (f)-(g), when l/Ueff intersects with xo> a true resonance 
occurs within the RPA. In (d)-(e), the 1/(7 in the Us vertex does not 
directly intersect to total x'o, but to its intraorbital components only 
and thus not shown. LI3j (c) The corresponding LDA FS in the folded 
zone is plotted at the same doping. The white to black background 
color gives the s*— symmetry. The two arrows point to the nesting 
channels which are responsible for the two modes in (a). The calcu- 
lation is done with SC gap Aq = A^/^ = 2Ap = 7meV. 

Double resonances: Fig. 1 shows our computed magnetic 
spectrum as a function of energy along the diagonal direction 
in the SC state for an representative doping of x = 0.20. The 
essential ingredient for the appearance of two resonances is 
having different amplitudes of gaps on different pieces of FSs. 
ARPES measures the average SC gaps which follows roughly 



for all dopings as Aq ~ ^-y/s ~ 2Ap, where a and /3 are the 
inner and outer hole pockets at F, respectively and 'y/S are the 
electron pockets at AI = (tt, tt) as shown in Fig. 1(c). 

The superconductivity is included in BCS formalism 
with s*— pairing symmetry as Ai(k) = Aj[cos(fc2:a) + 
cos {kya)]/2, where i is the band index. With these exper- 
imental inputs, we find a clear signature of two resonance 
modes: wjig^ ~ 0.75 x (2Aq) near the commensurate momen- 
tum Q — (tt, tt) and the second one at uj'^^^ ^ 0.5 x (2Ac) < 



w^gs at an incommensurate vector q ~ 0.78(7r, tt). This is the 
main result of this paper. 

According to the conventional view, apart from the ME ef- 
fect the magnetic structure in BCS Xo is entirely governed by 
the sign change of gaps at the 'hot-spots' and the energy con- 
servation formula for inelastic scattering of Bogoliubov quasi- 
paiticle on the FSlfT4ll 



(kp,q) = |ALi + iA;:' 



F+ql 



(2) 



At this locus on the (q, w)— phase space, x'o attains a logarith- 
mic jump, and due to the Kramers' Kronig relation Xo have 
a discontinuous peak as shown by the blue solid and green 
dashed lines in Figs, l(d-e). Therefore within the RPA, for 
a broad range of interaction U, a resonance is possible at the 
same locus. |9 | With these observations, apart from the true in- 
tensity of Xq, all its dispersion in the magnetic spectrum can 
be understood from the spanning vectors qs between various 
FS pieces which have opposite signs of the SC gap. 

At the doping considered in Fig. 1, all four pockets are 
present on the FS; see Fig. 1(c). For the case of s*— pairing, 
the gap below the magnetic BZ is always positive and above 
that, it is negative; i.e., a and /3 FS always have a positive 
sign while 7 and S bands have a negative sign at all dop- 
ings studied here, as illustrated by the white to black color 
background in Fig. 1(c). Thus only the inelastic scattering 
between a, /3 o 7, 5 is allowed, and the others including in- 
traband scattering are prohibited. The a— FS piece resides 
closer to the F— point than the f3 one, and •j/S lies close to 
M ~ (vr, tt) points. Therefore the corresponding magnetic 
vector Qi for the scattering channel a o 7, (5 is larger than 
Q2 for /3 o 7, (5. 

To clarify the physical origin of the two energy scales 
within the RPA, we take two representative q cuts for Xo, Xo 
and x"iRPA), Fig. 1(d) and (e). At Q = (tt, tt) in Fig. 1(d), 
Xo shows two characteristic peaks at uj^^^ and w^g^ coming 
from Qi and the residual Q2 scatterings respectively. The for- 
mer peak is stronger than the latter as the corresponding vector 
Qi spans up to Q but Q2 < (tt, tt). The situation is reversed at 
q = 0.78(7r, vr) where Qi > 0.78(7r, tt) but Q2 - 0.78(7r, vr), 
and thus the lower energy peak gains more intensity while the 
higher energy one goes negative, prohibiting the possibility of 
any resonance to occur within the RPA. 

In the RPA, the resonance peaks are set by the locus [Eq. 2] 
and intensity (ME in Eq 1) of the discontinuities of Xo- We use 
U = l.leV, a choice that is reasonable; moreover, the results 
are representative and remain qualitatively the same for dif- 
ferent U. The value of U is close to the value of 1 .2eV used 
in an earlier calculation to produce one resonance, which is 
closer to the first resonance we find here. 1 8 1 With this choice 
of U, the higher energy peak in Xo yields resonance peak at 
[Q,uj^^g], which is clearly separated from the second reso- 
nance at [0.78(7r, tt), iOresl- Away from these two sharp reso- 
nance peaks, the intensity is spread out over the large (q, oj) 
region, showing a remarkable downward dispersion as shown 
Fig. 1(a). The downward dispersion in pnictide can be con- 
trasted to the upward dispersion seen in chalcogenides. [15i 
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FIG. 2: (Color online) (a)-(c) Computed spectra of x" (q, i^) are plot- 
ted in color for underdoping in (a), optimal doping in (b) and over- 
doping in (c) of the hole-doped sample. The calculation is done here 
within in the minimal model described in the text. The magenta and 
yellow lines are the q and iu integrated value of x", respectively. 
The two energy scales in (b) match well with the available data at the 
same doping in Ref. 10 (d)-(f) The ARPES FSs corresponding to 
the INS spectra given in their upper panel at the same dopings||5j|. 



In the rest of the paper, we introduce a minimal model cal- 
culation of the INS spectra using ARPES or LDA FS which 
may still be valid in a multiband system even when the ME 
effects are neglected [M = 1 in Eq. 1]. We investigate the 
consequences of avoiding ME effect in Figs. 1(b) and l(f- 
g). x'o shows qualitative similarity between the calculation 
with and without including ME [compare Fig. l(d-e) with 
Fig. l(f-g)]: While both the peaks are still present, the lo- 
cations of their maxima are shifted slightly in both q and uj 
space. But the relative intensity of the peaks is subject to 
change, which can be compensated by tuning the effecting 
interaction ?7e//- Technically, setting the ME to be 1, we 
get all components of x tensor to be equal in Eq. 1. Ap- 
proximating the RPA interaction vertex Us by an effective 
scalar interaction Uef f < U, we can simplify the total RPA 

X = Exo(l - f/sXo)~^ w Xo(l - UeffXo)'^ (the sym- 
bol tilde over a quantity signifies that it is a tensor). We find 
that Ueff = O.leV can reasonably reproduce the results ob- 
tained with the ME effect including both resonances and the 
downward dispersion of the spectra. Of course this approxi- 
mation does not preclude us from analyzing the full problem 
but merely simplifies the analysis. 

In this approach, we compute the INS spectra with the RPA 
correction using the FS obtained either from ARPES when- 
ever available or from the LDA. The other input in our calcu- 
lation is the gap values which we obtain from various probes 
as listed in the corresponding figures below. An added benefit 
of this approach is that it allows us to compute the INS spec- 
trum entirely from experimental inputs, which is one of our 



FIG. 3: (Color online) Computed INS spectra plotted for two un- 
derdopings in (a,b) and optimal doping in (c) of the electron doped 
case, (d)-(f) The corresponding FSs are given at their lower pan- 
els at the same dopings obtained from the LDA band in (d) and (e) 
and from ARPES data] 16 1 in (f). The gap values in (d) and (e) are 
available from penetration depth measurement] 17], and optics llSI 
respectively at the same dopings where the same values for (c) are 
obtained from ARPES 1 16 1. In the LDA, FS pockets are known to be 
slightly larger than the experimental value due to the neglect of the 
spin density wave order. Hi 91 . 

goals. 

Doping dependence of two resonances: We compute the 
doping evolution of the two resonances by using experimental 
SC gaps and the ARPES FS. We note that while the doping 
dependence of the two resonances can also be obtained from 
the full BCS-RPA calculation (not shown), the energy and q 
values of the resonances are determined more accurately from 
the ARPES FS due to its inconsistency with LDA FS. 

Fig. 2 shows our computed spin resonance spectrum along 
the diagonal direction in the SC state as a function of dopings 
in hole-doped Fe-122 compounds. The separation between 
the two resonances in both energy and momenta is present for 
the doping range when both the hole pockets at F are present. 

The first mode, w.^es, occurs at Q = (tTjTt) at all dopings 
with a dispersion branch which is characteristic of the shape 
of the a— FSs. On the contrary, the doping dependence of 
the second mode Q2 is significantly large, deviating from a 
commensurate to incommensurate spectra as the contributing 
/3— FS grows in size with hole doping. We predict that with 
further underdoping, the second mode persists even when the 
first mode may disappear at a doping when the a— FS van- 
ishes (see electron doping below). Furthermore, we predict 
that with overdoping, when the electron pockets at M disap- 
pear, both these resonances will vanish. And, if the pairing 
is still present, the corresponding gap symmetry will change 
to the one for which the gap will change sign in the same FS 
pieces (such as d^2_y2 or dxy)', the neutron mode will be con- 
fined to the small q-space (see supplementary material). 
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For electron doping the q dependence of the second mode 
is opposite to the hole doping as the doping evolution of its 
FS topology is opposite to that of the hole-doped one, see 
Fig. 3. With increasing electron doping, the hole pockets at 
r gradually disappear while the electron pockets grow in size. 
The spatial separation between the two modes gradually de- 
creases as the distance between the a— and FSs decreases 
and both modes tend to be more commensurate in nature than 
that in the higher doping of the hole-doped side. This behav- 
ior is consistent with the reasons why the static commensurate 
spin-density wave order is more stable in electron-doping than 
in the hole doped one.||6l|7][T£j Finally, when the a band van- 
ishes from the FS, oj^g^ also vanishes; see Fig. 3(c). 

For electron doping, even when the first resonance van- 
ishes at x = 0.15, the second resonance is still present at 
ujres — 9.3meV at Q for the ARPES FS and the experimental 
gap amplitudes and agrees reasonably with the experimental 
value of 8.6meV at Q at a slightly larger doping of a; = 0.16 
[Ref. 12)]. And with further electron doping, the /3— FS also 
vanishes and the second resonance disappears. In this doping 
the superconductivity is seen to dis appear. Il20l 




Doping (x) 

FIG. 4: (Color online) Doping dependence of two resonance ener- 
gies and the SC gap amplitudes. The open symbols of various colors 
gives the SC gap values [red circles=Ac,, blue squares=A;3, green 
stars=A.^/i] as taken from experiments (see text). The filled sym- 
bols are the two resonance energies [magenta lower (cyan upper) 
triangles=aj^gs shown in Figs. 1,2. The shadings of same 

colors are a guide to the eyes. 

The two modes are not only q— resolved, they also appear 
at two energy scales as their underlying scattering states have 
different gap values. According to Eq. 2, each of the INS 
spectra is directly related to the SC gap values of partici- 
pating scattering states. For example, the resonance at Qi, 

= Aa|5£,l + A^i&\gC+Q, \ ~ 2A„5£^, by using the 

fact that Aa ~ A^^s and g^^ — —g^^f_^^- Therefore, in 
pnictide the leading or the higher energy resonance is always 
w^gj, cx 2Aa at all dopings (see Fig. 4). The proportionality 
constant is doping dependent, set by the value of the 'hot-spot' 
momenta. At optimal doping x = 0.4 of the hole doped case, 
we find that cj^^^ — 13.7meV and w^^^ — lOmeV by using 
the ARPES FS and gap values which agree remarkably well 
with the experimental value of 14 and 9.5meV, respectively at 
the same doping[10|. There is no one-to-one correspondence 



of the resonance energy with Tc as 2A/ ksTc is found to be 
doping dependent ||5J. 

This universal relationship between the resonance energy 
and SC gaps indicates that pnictides are in the weak- or 
intermediate-superconducting-coupling region where mean- 
field treatment of the SC state works well. Interestingly, the 
intensity of the resonance is found experimentally to follow 
the same T dependence of the SC gaps and vanishes at Tc in 
pnictide[ll| and chalcogenides|21 1. The universal relation- 
ship ujres fx (Ai + A2) (where Ai, and A2 are the SC 
gap amplitudes of the two orbitals connected by the 'hot- 
spot') is also found to be consistent with other families of 
superconductors. [ 1] 

Conclusion: In conclusion, we have shown that for a multi- 
band system with unconventional pairing with sign-changing 
gaps, one expects to find multiple spin resonances which are 
separated both in energy and momenta. The essential mech- 
anism for having two resonances in pnictide is the splitting 
between two hole pockets with different gap amplitudes. The 
difference in SC gaps and g— vectors of these two bands are 
the main reason why we see two resonances. Any calcula- 
tion invoking two bands will thus miss the second resonance. 
These results can be compared with INS experiments with the 
aim to observe second resonance in pnictides. 

The present computation does not have any free parameter 
as both the required FS information is obtained from either 
ARPES or the LDA and the SC gap values are taken from 
experiments. Our method of analysis reveals an intriguing re- 
lationship between the ARPES and INS measurements and 
one spectroscopic information can be extracted from the other 
one. This will give an important tool for any other multiband 
superconductors to study the pairing symmetry and/or the FS 
information which are still not clear 
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